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Introduction 


As a result of a study into the ability of modem Computational Fluid Dynamics (CFD) 
codes to adequately analyze aeroservoelastic interactions, it was concluded that there exists a 
strong need for well-documented experimental data which could be used to validate and/or 
calibrate these codes. Although extensive data is already available in the literature, much of it is 
unsuitable for use in this manner. In view of this, the Structural Dynamics Division of the 
NASA Langley Research Center (LaRC) has initiated an experimental program in aeroelasticity. 
The primary aim of the program, the Benchmark Models Program (BMP), is to provide 
experimental data to evaluate CFD codes for aeroelastic analysis. The program also aims to 
provide insight into the phenomenon of unsteady flow and to provide a data base for empirical 
design. (See refs. 1-4.) 

The BMP involves a series of aeroelastic models with different configurations to be 
tested in the Transonic Dynamics Tunnel (TDT) of NASA LaRC. One of these models is the 
Benchmark Active Controls Testing (BACT) model (ref. 1). It consists of a rigid wing of 
rectangular planform with a NACA 0012 profile and three control surfaces, namely, a trailing- 
edge control surface, a lower-surface spoiler, and an upper-surface spoiler. Figure 1 is a 
photograph of the BACT model, whilst figure 2 is an illustration of its planform. The model will 
be tested in the TDT using a flexible mount system which only allows plunging and/or pitching 
motion. Whilst having goals consistent with the BMP as a whole, a further important aim of the 
testing will be to validate analytical design tools for active control. 

One data set of interest to the program would be that which results from the model 
oscillating sinusoidally, in turn, in pure plunge and pure pitch. Of vital interest, then, would be 
the forces required to cause the model to execute the prescribed plunging or pitching motion. 
Clearly this is because the choice of the actuation system (for example, low force-capability gas 
thrusters or high force-capability hydraulic jacks) depends on the force requirements. It is with 
the quantification of these force requirements that this report is primarily concerned. 

It will be noticed that if the forces required to move the model could be determined 
exactly, it must mean that the aerodynamic forces acting on the model are known exactly, and 
hence the need would not exist to perform a test to measure these aerodynamic forces. (A major 
aim of the BMP is to determine the forces acting on wind tunnel models.) It should therefore be 
emphasized that the purpose of this report is to determine approximately the forces required to be 
applied to the model, so that the feasibility of this part of the BACT model test may be assessed. 
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Symbols 


NOTE:(l) The term “wing" is meant to include the effect of the (fixed) trailing-edge control 
surface. For the purposes of this analysis, the spoilers are non-existent. 

(2) The elastic axis, center of mass and mid-chord of the wing are coincident. 

a speed of sound in tunnel medium 

« « 22/«77 

b semi-chord 

b (A;+A2)/n/7 

c ^ih/(nj]n22 -nnmi ) 

dA distance of Fa from mid-chord of wing 

ds distance of Fb from mid-chord of wing 

d c distance between trailing-edge control surface hinge line and center of mass of trailing- 

edge control surface 

dh distance of trailing-edge control surface hinge line from mid-chord of wing 

/ frequency 

Fa external point force acting on model at distance d\ forward of mid-chord of wing 
(positive up) 

Fb external point force acting on model at distance ds aft of mid-chord of wing (positive up) 
fi wind-off natural frequency 

f(x) aerodynamic force per unit chord length 

gi I th modal damping ratio 

H hinge moment about trailing-edge control surface hinge line (positive nose-down) 

i V=T 

/ moment of inertia of wing about mid-chord of wing 
I identity matrix 

/<£ moment of inertia of trailing-edge control surface about its center of mass 

k plunge stiffness of pitch and plunge mount 

k reduced frequency, cob/V 

K pitch stiffness of pitch and plunge mount 

K generalized stiffness matrix 

ki I th reduced frequency 

kij ij d* term of K 

L lift force acting through mid-chord of wing (positive up) 
m mass of wing 
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M 

M 

mij 

n ti 

P 

q 

Q 

QD 

<W 

Wjl 

WjR 

Sij 

t 

T 

U 

V 

V i 

vg 

X 

y 

Y 


Mach number 
generalized mass matrix 
ij d 1 term of M 

mass of trailing-edge control surface 
ij element of M _1 

pitching moment about mid-chord of wing (positive nose-down) 

generalized displacement vector, given by z=[ y 6 £ ] T 

matrix of generalized aerodynamic force coefficients 
1 2 

dynamic pressure, ^pV 

given by QDiQijp^Qiji)’ that is, the ij^ term of <?dQ 

imaginary part of y* term of Q 

real part of y * term of Q 

contribution of y th mode to i th sensor reading 

time 

kinetic energy of system 
potential energy of system 
free-stream velocity 
gust velocity (positive up) 

non-dimensional gust velocity, given by vg=v g /V (positive up) 
chordwise distance from mid-chord of wing (positive aft) 
vertical displacement of mid-chord of wing from equilibrium (positive down) 
amplitude of plunging motion 


virtual work done by aerodynamic forces 
6 pitch angle of wing from equilibrium (positive nose-up) 

& amplitude of pitching motion 

X eigenvalues of M _ 1 K 

A,- i * non-zero eigenvalue of M" 1 K 

% deflection of trailing-edge control surface relative to wing (positive trailing-edge down) 

p density of wind-tunnel medium 

CO circular frequency 

C0y circular frequency of plunging motion 

t oq circular frequency of pitching motion 
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Equations of Motion 


The equations of motion will be derived to take into account the effect of trailing edge 
control surface motion and vertical gusts. This is because the numerical model from which the 
aerodynamic forces are extracted include the forces acting on the wing due to these two 
phenomena. The spoilers, for the purposes of this analysis, are non-existent. 

Prior to any testing, the model will be balanced so that the system center of mass 
coincides with the elastic axis of the mount system, which also corresponds to the mid-chord 
position. This decouples the plunging and pitching modes, provided the trailing-edge control 
surface is fixed. (Compare with ref. 4.) 

For the purposes of this analysis, the effect of gravity is ignored. Although this effect is 
of great importance as far as the position of the trailing-edge control surface is concerned, 
neglecting it does not affect the results as the trailing-edge control surface is later fixed. 

The equations of motion may be formulated by making use of Lagrange’s equations (ref. 
5). With reference to figure 3, the system kinetic energy and potential energy are respectively 
given by 


T = +\ m t;[y+( d h+d c )e+ d c $f + [e + $ f 
-\ m l;[y+{ d h +d c)6f and 

U = -lcy 2 +-Kd 2 . 

2 2 


0) 


( 2 ) 


Using Equations (1) and (2) in conjunction with Lagrange’s equations gives, as the in vacuo 
equations of motion, 


Mq + Kq = 0, where 


( 3 ) 


M = 


m 


0 


0 

I 


m S d c 


m s( d h+ d c) d c 


+n 


j % i d h + d c ) d c , 2 . i 

m $ d c + / e mgdc+Ig 


I 
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y 

e . 

4 . 

It will be noted that the effect of structural damping has been neglected. This is justified 
by the fact that any structural damping present will be small. (Table 1 of ref. 4 quotes an 
experimentally measured structural damping figure of 0.0024 for both the plunge and pitch 
modes for a NACA 0012 aerofoil with the same mounting system.) 

Now let the system be acted upon by external non-conservative forces. These forces 
consist of (1) the aerodynamic forces, and (2) any other external non-conservative forces. The 
aerodynamic forces are represented in figure 4. The virtual work done by these forces is given 
by 



k 

0 

O' 


K = 

0 

K 

0 

and q = 


_0 

0 

0_ 



8W = -j* b f(x)dx(8y + x89 ) - f(x)dx(x - d h )8£ 

= Sy^ b f(x)dx - 89 \ b _ b xf{x)dx- 8$ (* - d h )f(x)dx. 


That is. 


8W = (-L)8y + (-P)86 + (-H)8Z, where (4) 

L = t b f( x ) dx ’ 

P = \ b xf(x)dx, and 
J-b 

H = j b ^(x-d h )f{x)dx, 


are the aerodynamic lift, pitching moment (about the mid-chord of the wing) and hinge moment 
(about the trailing-edge control surface hinge line) respectively. Hence, the equations of motion, 
including the aerodynamic forces, become 

Mq + Kq + [L P H] J = 0. (5) 

Now, for constant amplitude sinusoidal motion , the vector [L P H ] T , the generalized 
aerodynamic force vector, may be represented by 
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(6) 


[L P H} 1 = q D <l\y e t; v g ] T , 

''g =v g /V bein g ^e vertical gust velocity, non-dimensionalized by the free-stream velocity, 

v > Q a matrix of generalized aerodynamic force coefficients, and qo the dynamic pressure (ref. 
6). (The contribution due to gusts is later excluded from consideration.) 

Now, if, as shown in figure 5, the external point forces F A and F B also act on the 
structure, the complete equations of motion, including the effect of the external non-conservative 
forces, are 


Mq + Kq + q D Q 



+ 


1 1 

~ d A d B 
0 0 


Fa 

f b 


= 0 . 


(7) 


Determination of External Forces 


Having derived the equations of morion in the form of equation (7), it is then possible to 

perform an inversion to determine the point forces, F A and F B , in terms of the generalized 
coordinates.* 


Let the ij lh term of the generalized mass, stiffness, and aerodynamic force coefficient 
matrices be given by mq, kq and qq R + iqq t respectively. Setting £ = v g = 0 (that is, fixing the 


trailing edge control surface and removing the effect of gusts) with 


y = Tsin (Oyt and 
6 = Ssino) B t 


gives 


[~™1 1 co y 2 + ki i + q D 1r + iq x l{ jj Y sin (o y t 
+ QD (qi2 r + *<712/ ) ® sin 0) e t + F a + F B = 0, and 


In this sense, the problem described herein is a linear Inverse Problem, since the unknown input to a known 
system with a known response to this input is sought. In fact, the problem falls into the category of the 
Reconstruction Problem (ref. 7), or, more precisely, the Force Identification problem. 
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( 10 ) 


<7d(< 7 21* +“721 / ) } ' sinG y 

+[ _m 22 COq 2 + k 2 2 + qo(q 22 R + “722/ )] ® sin 03 6 l ~ p A d A + F B d B = 0. 

as the equations of motion, in which plunging and pitching motion occur with amplitudes of Y 
and 0, and with circular frequencies of a>y and cog respectively (t = time).* Solving equations (9) 
and (10) simultaneously for F A and Fb gives 


p d B 

r a — 


dA + d. 


B 


2 1 1 
— m\\C0y + &n+<7n <721 I Y sin 6)yt 

ds ) 


1 


1 


1 


-J-m H^e 2 + <7l2 -~r k 22 ~~qi2 
d B d B d B 


©sinty^r 


( 11 ) 


, and 


f b = 


- Aa 

dA+d B 


LV 


-m u (O y 2 + k u +<7n +— <721 J 


Y sin CQyt 


( 


1 2 1 , 1 U . 

m 2 26>0 +<?i2 + — ^22 + 3~ ^22 ©smOgf 
d A d A d A J 


( 12 ) 


where qij = qoiqijR + iqijj), namely, the i/ th term of the matrix of generalized aerodynamic force 
coefficients Q, multiplied by the dynamic pressure qo- The point forces, F A and F B , required for 
the model to execute pure plunging and pure pitching motion may then be easily determined as 
follows. 


Pure Plunge 


If the model executes pure plunging motion, the motion of the model as a function of 
time may be described by the equations 


y = ysinty^r and 

0 = 0. 


(13) 


Thus, 0 = 0 is substituted into equations (11) and (12) to yield 


* The amplitudes may be complex to represent arbitrary phase differences between the plunging and pitching 
motions. 


7 



f a = 


-d B y 

d A +d B 


-m n Q) y 2 + k\\+q D q UR - ~~ q D q 2iR j 


+ i 


qD q\\] -~^qDq2\i 


sinc^r , and 


( 14 ) 


f b = 


-d A Y 


dA +d B 
( 


( 2 1 N 

m i l My + *i l + <7z> <7i t* + — <7 d q 2 i R 

A 


+ Z 


1 


^£>< 711 / +—qDq2\j 


sin (Oyt 


(15) 


as the forces required for the model to execute pure sinusoidal plunging motion. 

Pure Pitch 

Similarly, if the model executes pure pitching motion, the motion of the model may be 
described by 


y = 0 and 
6 = ©sinoigr. 


(16) 


Hence, Y = 0 is substituted into equations (1 1) and (12) to give 


( 


1 


1 


—m 22 0) e + q D q ]lR - —k 22 - —q D q22 


+i 


\qDq\2j ~—qoq22, 


d B d B 
sin cogt, and 


(17) 


r -d A e 
F B = — — 


d A 

( 


m 22 0 )Q 2 + <?£> q\2 R + —*22 + ~J~qDq22 R 


+ l 


1 


qoq 12 / +— qoq22j 


sin (Ogt 


( 18 ) 


as the forces required for the model to execute pure sinusoidal pitching motion. 
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Numerical Results 


With the aid of equations (14), (15), (17) and (18), the magnitudes (per unit of plunge or 
pitch amplitude) and phases (relative to the motion) of the forces required to cause the model to 
execute each of the two types of motion may be determined. The numerical values of the various 
terms in these equations may be determined as described below. 

The Interaction of Structures, Aerodynamics, and Controls (ISAC) computer program 
(ref. 8) was previously used to perform an aeroservoelastic analysis of the BACT model without 
spoilers. It is from the results of this analysis that the majority of the numerical values required 
are obtained. 

It should be noted that the ISAC numerical model is three-dimensional as far as the 
aerodynamics are concerned (that is, there is a spanwise and chordwise distribution of 
aerodynamic forces). This three-dimensional model is reduced to a two-dimensional one by 
making use of the generalized aerodynamic force coefficient matrix, Q, which contains the 
forces (normalized by the dynamic pressure qo) acting on the two-dimensional model described 
herein due to changes in the coordinates y, 6, £ and vg used to quantify this two-dimensional 
model. 

Appendix A contains a listing of the ISAC output file DYN_TAPE5.DAT for a Mach 
number, M, of 0.78. The format of this file is illustrated in appendix B. As may be seen, the 
generalized stiffness matrix is not explicitly given. This matrix may be determined from the 
generalized mass matrix and the natural frequencies by the following method. 

Let A represent the eigenvalues of M**K, and hence the square of the natural frequencies 
(in rad s* 1 ) of the undamped system described by equation (8). Then, detfM^K - A I) = 0, or 


det 


«i \k - A 

»2l* 

« 31 * 


n n K 

«22^ — A 

n 22 K 



where nij is the ij& element of M' 1 , which in turn leads to 

A^A 2 ~ {n\\k + n 22 ^)X + (nn«22 _n 12 n 2l)^] = 0- 
If 0, Aj and A 2 are the eigenvalues. 


(19) 


( 20 ) 


A 2 — {/i\\k + n22^)X + (nii«22 ~ n l2 n 2l)kK = A 2 — (Aj + A2)A + A 1 A 2 . (21) 
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Equating coefficients and solving for k and K gives 


k = 



and K~ — , 
k 


( 22 ) 


where a — «22/ n ll> ^ — a °d c ~ ^-1^2 ^ i n l\ n 22~ n l2 n 2l)- As may be seen from 

appendix A, the ISAC model has plunge and pitch natural frequencies of 3.40 Hz and 5.16 Hz 
respectively. Thus, set = (3.40 x 2 iCp- and A 2 = (5.16 x 2 lip- (rad^ s'^), and then compute M'l 
to obtain 


k = 21 6.459, 498.555 lb in- 1 and 
K = 33008.2, 14331.3 in-lb rad-1. 

An examination of the eigenvectors of M-lK for each solution pair (k, K) (see fig. 6 for a 
diagrammatic representation of these eigenvectors, or, alternatively, the mode shapes of the 
system) reveals that the solution pair (k = 216.459, K = 33008.2) results in plunge occurring 
before pitch, whilst the solution pair (k = 498.555, K = 14331.3) results in pitch occurring before 

plunge. Since the plunge frequency (3.40 Hz) actually occurs before the pitch frequency (5 16 
Hz), 


£ = 216.459 lb in'l and K = 33008.2 in-lb rad’ 1 

should be chosen as the applicable stiffness values. 

The generalized aerodynamic force coefficients, qij R + iqij , , may also be obtained from 

the ISAC output file of appendix A. Contained in appendix B is an explanation of the manner in 
which these are listed. 

Finally, numerical values are required for d/^ and dg . For the purpose of this analysis, 
these will be taken as 


dA = ds = 7 in. 


Using the above-mentioned numerical values, the FORTRAN program PPRFCE in appendix C 
was written to determine, from equations (14), (15), (17) and (18), the magnitudes (per unit of 
plunge or pitch amplitude) and phases (relative to the motion) of the point forces. Fa and Fb, 
required to constrain the model to execute each of the two types of motion (at the frequencies 
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corresponding to the reduced frequencies for which the generalized aerodynamic force 
coefficients are determined by ISAC). The tunnel medium is Freon, in which the speed of sound, 
a ~ 500 ft s-l. (This is for a 95% Freon 12/air mixture at a temperature of 530* R and a pressure 
of 2000 psf.) Using this value for a leads, with a model semi-chord, b , of 8 in, to the relationship 
between frequency, reduced frequency and Mach number illustrated in figure 7.* 

Figures 8-11 contain the plotted results of four separate runs of PPRFCE for the 
following test conditions: 

(1) Wind-off 

(2) M = 0.78, q = 100 psf (below flutter) 

(3) M = 0.78, q = 146 psf (near* flutter) 

(4) M = 0.78, q = 200 psf (above flutter) 


Discussion 


The purpose of this analysis is to approximately determine the forces required to move 
the BACT model in pure plunge and pure pitch. The graphical data contained in figures 8-1 1 
presents this information for various test conditions. 

As expected, the wind-off behavior of the system is that of a classical undamped two 
degree-of-ffeedom system. As may be seen from figures 8(a) and 8(b), no force is theoretically 
required at either the plunge (3.40 Hz) or pitch (5.16 Hz) natural frequencies, to produce, 
respectively, pure plunging or pure pitching motion, with the forcing function shifting from 
being in phase with the motion to being out of phase with the motion as the frequency of 
excitation is increased through each natural frequency. (The lack of frequency resolution in the 
ISAC data used prevents this phenomenon from being exactly represented graphically in figs. 8- 
11. Also, note that because of the directions defined to be positive for Fa, Fg,y and 6, a 180° 
phase angle between Fa and F g and the plunge coordinate y, for instance, means thatF^ and Fg 
are actually in phase with y.) 


From this figure, a connection may be drawn between the forces required to move the BACT model (which, as 
will be seen later, can be highly frequency dependent) and the reduced frequency achievable with these forces. A 
great deal of emphasis is placed on this reduced frequency parameter for the BACT model test, 
t The previously-determined flutter dynamic pressure at M = 0.78 is 146.135 psf. This was obtained by using the 
STABCAR module of ISAC. 
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For the wind-on test cases, it will be seen that the effect of aerodynamics on the forces 
required to produce plunging motion is small (compare fig. 8(a) to figs. 9(a), 10(a) and 11(a)). 
The effect of aerodynamics is, as might be expected, greater when the model executes a pitching 
motion (compare fig. 8(b) to figs. 9(b), 10(b) and 11(b)). This may especially be said to be so 
considering that the plunge forces are shown per inch of plunge amplitude, while the pitch forces 
are shown per degree of pitch amplitude. Because the center of pressure is closer to the point of 
application of Fa than it is to that of Fb, Fa is affected to a greater extent than Fg. As my be 
seen from figures 9(b), 10(b) and 11(b), the minima of the IF^I vs. /and IF^I vs./ curves for 
pitching motion now occur at different points, with the result that there is no one frequency 
where both the forces required to achieve pure pitching motion are theoretically zero. In fact, for 
the test conditions considered, the quantity [(IF^I, IFflOmax] m j n is of the order of 20-30 lb per 
degree of pitch amplitude, and occurs when the curves intersect. In other words, the force 
actuators which are used to move the BACT model in pure pitch must be capable of producing 
20-30 lb per degree of pitch amplitude before any pure pitching motion is possible. As with 
pure plunging motion, this minimum force requirement (which is theoretically near zero for pure 
plunging) occurs at a particular frequency, and increases as the amplitude of motion is increased, 
the frequency of motion is reduced, or, most severely, the frequency of motion is increased. 

It may therefore be seen from figures 8-1 1 that, for the test conditions considered, only 
plunging motion may be achieved with negligible force. This is possible, wind-off and wind-on, 
only near the plunge natural frequency of 3.40 Hz. (For M = 0.78, this, according to fig. 7, 
results in a reduced frequency, k = 0.035.) As far as pure pitching is concerned, the smallest 
force, per degree of pitch amplitude, that will produce the motion is of the order of 20-30 lb. 
This occurs near/= 4 Hz (k = 0.04). Excitation at higher frequencies, in either plunge or pitch, 
requires significantly larger forces. This is attributed to the quadratic relationship between the 
inertial loads and the frequency of motion. 


Concluding Remarks 

Presented above is an approximate analytical determination of the forces required to 
move the BACT model in either pure plunge or pure pitch. The analysis makes use of results 
previously obtained from ISAC to obtain the majority of the required numerical values. 

It was seen that, for the test conditions considered, only plunging motion may be 
achieved with negligible force. It was pointed out that this is possible, wind-off and wind-on, 
only near the plunge natural frequency. As far as pure pitching is concerned, the smallest force, 
per degree of pitch amplitude , that will produce the motion is of the order of 20-30 lb. This 
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occurs near a frequency slightly less than the pitch natural frequency. It was also pointed out that 
excitation at other frequencies, in either plunge or pitch, requires larger forces. This increase in 
the force requirement is especially significant at higher frequencies. 
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Appendix A— ISAC DYN_TAPE5.DAT File for M = 0.78 


ASEB2> sh def 

525DUC0 : [DCRUZ.PAPA RIG.MATX. ISAC_DATA] 
ASEB2> ty dyn_tap«5 m7s.dat 

0.000000 O.OOIOOTJ 0.005000 0.010000 

0.100000 0.130000 0.160000 0.200000 

0. OOOOOOOE+OO 0. 0000000E+00 
0. OOOOOOOE+OO 0 .OOOOOOOE+OO 
0 . 0000000E+00 0. 0000000E+00 

0.2360745E+04 0 .OOOOOOOE+OO 

-0. 1060751E+05 0 .0000O00E + 00 

0 . 4733114E+02 0. 0000000E + 00 

0.4571392E+03 0 . 0000000E+00 

0 . 3123492E+03 0 . 0000000E+00 

0 . 1047493E+03 0 . OOOOOOOE+OO 

0 . 2360745E+04 0 . OOOOOOOE+OO 

-0.1060751E+05 0 .OOOOOOOE+OO 

0.4733114E+02 0 . OOOOOOOE+OO 

0 . 2086285E-03 0 .2950917E+00 

-0.2832395E-02 -0 . 1325931E+01 
-0.2807885E-04 0 . 5916387E-02 


0.020000 0.040000 0.060000 0.080000 


0.2360736E+04 
-0 . 1060746E+05 
0 . 4733095E+02 
0.4571370E+03 
0.3123663E+03 
0.1047494E+03 
0.2360725E+04 
-0 . 1060739E+05 


’9 . 3429534E+00 
0 . 2604929E+02 
0.4179317E+00 
*0 . 6140946E+00 
0 . 4 82O09OE+0 1 
0 . 1923599E+00 
■0 . 5355827E+01 
0.2987640E+02 


0.4733103E+02 -0 .1600142E-01 
0 . 5209487E-02 0 . 1475285E+01 

-0.7077803E-01 -0 . 6628576E+01 
”0 . 7020520E-03 0 . 2958 149E-01 

0 . 2360529E+04 -0 . 1705897E+01 
-0. 1060611E+05 0.1302027E+03 

0.4732767E+02 0 . 2089B02E+01 

0.4570744E+03 -0 . 3068 403E+01 
0 . 3127788E+03 0 . 2409324E+02 

0. 1047496E+03 0. 9618196E+00 

0 . 2360265E+04 -0 . 2676594E+02 
-0 . 1060470E+05 0 .1493114E+03 

0.4733041E+02 -0.7988867E-01 
0 . 2076123E-01 0 . 2949507E+0 1 

-0.2827187E+00 -0.1325045E+02 
-0.2809189E-02 0 . 5916107E-01 

0.2359888E+04 -0 . 3357057E+01 
-0. 1060193E+05 f. 2601355E+03 
0.4731824E+02 0 . 4180474E+01 

0.456B808E+03 -0 . 6124053E+01 
0 .3X40604E+03 0 . 4811731E+02 

0.1047505E+03 0 . 1923766E+01 

0 . 2358841E+04 -0 . 534504 6E + 02 
-0.1059630E+05 0 .2981854E+03 

0. 4732938E+02 -0 . 1590634E+00 
0 . 8186422E-01 0 . 58 90697E+01 

-0.1124816E+01 -0 . 2644846E+02 
-0 . 1125185E-01 0.1183044E+00 

0 . 2357406E+04 -0 . 6294 385E+0 1 
-0.1Q58S62E+05 0 . 5182003E+03 

0.4727958E+02 0 . 8367682E+01 

0 . 4561239E+03 -0 .1215026E+02 
0 . 3190988E+03 0 . 9570166E+02 


0 . 1047543E+03 
0 . 235324 4E+04 
-0 . 1056321E+05 
0.4732447E+02 
0.3109684E+00 


0 . 364B495E+01 
-0.1062712E+03 
0 . 5929800E+03 
“0 .3126648E+00 
0.1172026E+02 


“0 . 4413990E+01 
-0.4521089E-01 
0.2348475E+04 
-0.1052523E+05 
0.4714023E+02 
0 . 4533213E+03 
0.3381290E+03 
0.1047727E+03 
0 . 2332174E+04 
-0.1043750E+05 
0.4732190E+02 
0.6501896E+00 
-0 . 9671860E+01 
-0 . 1023003E+00 
0.2336057E+04 
-0 . 1043665E+05 
0 . 4694393E+02 
0 . 4492054E+03 
0 . 3670452E+03 
0 . 1048106E+03 
0.2300429E+04 
-0 . 1024525E+05 
0.4735741E+02 
0 . 1058729E+01 
-0.1667593E+02 
-0 . 1829174E+00 
0 . 2321967E+04 
-0.1032872E+05 
0.4671730E+02 
0 . 4441766E+03 
0.4037290E+03 
0.1048725E+03 
0.2260737E+04 
-0.1000096E+05 
0.4745755E+02 
0.149B740E+01 
-0 . 2521163E+02 
-0 . 2873B24E+00 
0.2307489E+04 
-0.1020785E+05 
0.4647630E+02 
0.4385271E+03 
0 . 44 65220E+03 
0 . 1049593E+03 
0 . 2215281E+04 
-0.9716916E+04 
0.4763315E+02 
0 . 214 7011E+01 
-0 . 404 6722E+02 
-0 . 4893222E+00 
0.2287047E+04 
-0 . 1001191E+05 
0 . 4611634E+02 
0. 4293600E+03 
0 . 5192844E+03 
0 . 1051388E+03 
0 . 21 40059E+04 
-0 . 9238680E+04 
0 . 4 804606E+02 
0 . 2705663E+01 
-0 . 58 18721E+02 
-0 . 7460511E+00 
0 . 22 69994E+04 
-0.9806817E+04 
0.4579297E+02 
0 . 4198221E+03 


-0.5250377E+02 
0 .2365471E+00 
-0 .9684481E+01 
0 . 1022035E+04 
0.1678198E+02 
-0.2362152E+02 
0.1876428E+03 
0.7703010E+01 
-0 . 2080625E+03 
0 . 1161 560E+04 
-0 . 5902 80 1E+00 
0.1744 777E+02 
-O.7707333E+O2 
0 . 354 938 6E+00 
-0 .081O893E+O1 
0 . 1504 642E+04 
0.2526579E+02 
-0.3408465E+02 
0.2737661E+03 
0 . 11 56420E+02 
-0 . 3028361E+03 
0.1690877E+04 
-0.8263804E+00 
0 . 2305262E+02 
-0.1023674E+03 
0.4739526E+00 
-0 . 34 95844E+0 1 
0 . 1964 955E+04 
0 . 3382457E+02 
-0.4346958E+02 
0 . 353294 5E+03 
0 . 1542908E+02 
-0.3896496E+03 
0 .2174084E+04 
-0 . 1035657E+0 1 
0 . 2852 791E+02 
-0.1258700E+03 
0 . 5941658E+00 
0.59821B9E+01 
0.2404241E+04 
0 . 4246000E+02 
-0 .51802 4 7E+02 
0.4259678E+03 
0 . 1929399E+02 
-0 . 4683183E+03 
0 . 2608 437E+04 
-0.1237728E+01 
0.3650571E+02 
-0 . 1591340E+03 
0 .7781810E+00 
0.2702299E+02 
0.3028612E+04 
0.5555920E+02 
-0.6247355E+02 
0.5221207E+03 
0 . 2508398E+02 
-0.5715457E+03 
0.3167B09E+04 
-0.1570158E+01 
0 . 4423395E+02 
“0 . 1899548E+03 
0.9687909E+00 
0 . 5468544E+02 
0.3619389E+04 
0 . 6883688E+02 
-0 . 7121 078E+O2 


0.5992074E+03 
0 . 1053749E+03 
0 . 2060313E+04 
-0 .8722548E+04 
0 . 4862977E+02 
0.3214034E+01 
-0 . 84 947 65E+02 
-0. 1174168E+01 
0.2254662E+04 
-0.9526773E+04 
0 . 4548124E+02 
0.4070796E+03 
0 . 7124526E+03 
0.1057707E+03 
0 . 1952159E+04 
-0 . 801 1389E+04 
0.4962676E+02 
4 . 7583333E-01 
0. OOOOOOOE+OO 
3 . 6212015E-03 
3.40 
5.16 
0 . 

0 . 0 . 0 . 

-1. -5. 0. 

-1. 7.15 0. 
ASEB2> 


0 . 6033608E+03 
0 .3085570E+02 
-0 . 6504796E+O3 
0 .3621006E+04 
“0 . 1990590E+01 
0 . 5423073E+O2 
-0.2272798E+03 
0.1237594E+01 
0.9873293E+02 
0 . 4 370152E+04 
0 . 8602436E+O2 
-0 . 8045405E+02 
0 . 6900458E+03 
0 . 3850715E+02 
-0.7526026E+03 
0 . 407501 4E+04 
-0 . 27 6924 8E+01 
0 . OOOQQOE+OO 3 
3 . 1 43999E+0 1 1 

1 . 790215E-02 8 


. 6212015E-03 

.790215E-02 

.575828E-03 
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Appendix B— Format of ISAC DYN_TAPE5.DAT Files 



£l 82 83 } Modal damping ratios 


Contribution of mode to i* sensor reading: 

(sij,j = 1,2,3), i‘ = 1.2 


(Units: slinch, in, s, rad) 
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Appendix C — FORTRAN Source Code 


ASEB2> sh def 

$25DUC0 : {DCRUZ . PAPA_RIG . MATX) 

ASEB2> ty pprfce.for 

PROGRAM PPRTCE 

THIS PROGRAM COMPUTES THE FORCES REQUIRED TO MOVE 
THE BENCHMARK ACTIVE CONTROLS MODEL IN EITHER 
PURE PLUNGE OR PURE PITCH. THE EFFECT OF 
AERODYNAMICS IS INCLUDED BY USING DATA GENERATED 
FROM ISAC RUNS. 


IMPLICIT NONE 

REAL'S PI,M,A,RFRQ(2Q),QR(20,20,20) , QI {20, 20, 20) ,GM<20,20), 
♦NFRQ(20) ,D,B,MU,M22,K11,K22,0M, FAPLR, FAPLI , F8PLR, FBPLI , 
+FAPIR, FAPII, FBPIR, FBP II , FAPL, FBPL, FAPI, FBPI, FREQ, FM(20> , 

♦F AFBP L ( 2 0 , 2 ) , FAFBPI (20, 2) , DUMMY, DYN, PHSEPL(20, 2) , PHSEPI {20, 2) 
INTEGER* 4 NX, I , NF, NC, NM, NG, J, K 

OPEN ( UNIT- 1, FILE- 'S2SDUC0: [ DCRUZ . PAP A_RIG .MATX . ISAC_DATA) 
+DYN_TAPE5_M78 .DAT’ , STATUS-' OLD ' , READONLY) 

OPEN (UNIT-2, FILE-' FAFBM7 8 QIOO.DAT* , STATUS- 1 UNKNOWN * ) 

C 

C READ IN ISAC OUTPUT: 

C 

PI-DACOS (-1.0D0) 

M-0.78D0 

A-500 . 0D0* 12 . 0D0 
NX- 12 

READU,*) (RFRQ(I) , 1-1, NX) 

NF— 2 
NC-I 

W-NF+NC 

NG-l 

DYN— 100 . ODO/14 4 . 0D0 
DO 30 K-l , NX 

DO 20 J-l , NM+NG 
DO 10 I— 1 , NM 

REMILM QR (I, J, X) , QI (I, J, X) 

QR ( I , J , K ) -QR ( I , J, K) *DYN 
QI (I, J, K) -QI { I, J, X) *DYN 
10 CONTINUE 

20 CONTINUE 
30 CONTINUE 

READ(1,*) ( < GM { I , J ) , J— 1 , NM ) , I - 1 , NM) 

READU,*) (NFRQ(I) , I-1,NM) 

C 

C CALCULATE FA AND FB FOR PURE PLUNGE AND 

C PURE PITCH: 

C 

D-7.0D0 
B-8.0D0 
M11-G1 (1,1) 

M22-GM (2,2) 

KU-216. 459D0 
K22-33008.2D0 
DO 90 I— 1 , NX 

OM-RFRQ(I) *M*A/B 

FAPLR— { -Mil *OM*OM+Kll*QRU, 1, 1) -QR (2, 1 , I) /D) /2 . 0D0 

FAPLI— (QIU, 1, I) -QI (2, 1, 1)/D)/2.0D0 

FBPLR-- {-Ml 1 *OM*OM+Kl 1+QR (1,1,1) +QR (2,l,I)/D)/2. 0DQ 

FBPLI— (QIU, 1, I) +QI (2, 1 , 1 ) /D) / 2 . 0D0 

FAPIR— (M22*OM*OM/D-K22/D+QR(l, 2, I)-QR(2, 2, I) /D) /2 .0D0 

FAPII— (QI (1,2, I) -QI(2,2, I)/D)/2.0D0 

FBPIR— (-M22*OM*OM/D+K22/D+QR U , 2, 1) +QR (2, 2, I) /D) /2.0D0 
FBP 1 1— (QI (1,2, I) +QI (2,2,I)/D)/2. 0D0 
FAPL—DSQRT (FAPLR*FAPLR+FAPLI ‘FAPLI ) 

FBPL-DSQRT (FBPLR*FBPLR+f BPLI ‘FBPLI ) 


TAPI-DSQRT <FAPIR*FAPIR+FAPII*FAPII) ‘PI/180.0D0 
rBPI— OSQRT (FBPIR*FBPIR+FBPII *FBPII) *PI/180.0DO 
FREQ-OM/ (2 .0D0*PI ) 

FM (I ) -FREQ 

FAFBPL (1,1) -FAPL 

FAFBPL (I, 2) -FBPL 

FAFBPI (I, 1)-FAPI 

FAFBPI (I, 2) -FBPI 

PHSEPL (I, 1 ) -DAT AN2D (FAPLI, FAPLR) 

PHSEPL (1,2) -DATAN2D (FBPLI, FBPLR) 

PHSEPI (I, D-DATAN2D (FAPII, FAPIR) 

PHSEPI (1,2) -DATAN2D (FBP II, FBPIR) 

90 CONTINUE 

CALL MATSAV ( 2 , ’MTXFM' , 20, NX, 1, 0, FM, DUMMY, * (1P3E25.17) *) 

CALL MATSAV (2, 'MTXFAFBPL * , 20, NX, 2 , 0, FAFBPL, DUMMY, ' (1P3E25.17) l ) 
CALL MATSAV (2, 'KTXFAFBPI 20, NK, 2, 0, FAFBPI , DUMMY, ' (1P3E25.17) ’) 
CALL MATSAV (2, ’MTXPHSEPL* , 20, NK, 2 , 0, PHSEPL, DUMMY, ’ (1P3E25.17) *) 
CALL MATSAV (2, • MTXPHSEP I * , 20, NK, 2 , 0 , PHSEPI , DUMMY, * (1P3E25.17) *) 
STOP 
END 

A5EB2> 
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Figure 1. BACT model with deflected aileron and deployed upper-surface spoiler. 



Figure 2. BACT model planform (reproduced from ref. 1, fig. 3). 
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Equilibrium 
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(a) k = 216.459 lb tar 1 , K= 33008.2 in-lb rad' 1 



Equilibrium 0 Hz 3.40 Hz 5.16 Hz 

(b) k = 498.555 lb in* 1 , K= 14331.3 in-lb rad" 1 


l igure 6. Mode shapes for different spring stiffnesses. 
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Figure 9(b). Pitch force characteristics (Af = 0.78, q = 100 psf). 


24 







Magnitude of 
pitch force per 
degree of pitch 
amplitude g 0C 

(lb deg 1 ) 


Phase of pitch 
force relative to 


Frequency (Hz) 














REPORT DOCUMENTATION PAGE 


Form Approved 
OMB No 0704-01 88 


f'^Di : reoorv^c tsuroe^ re mis rcmecrion •> *st r-'amc 

gamerira ana Maintaining the aata neeoea ana corroietm; ana r«.<® 
collection information, <nciuC!ng suggestions + or reoucing mis ou'o 
Davis highway. Suite 12C4 Arlington, „ « 22202-43C2 anc K f'® 0*^' 


1 AGENCY USE ONLY ( Leave blank) 2. REPORT DATE 

July 1993 


4. TITLE AND SUBTITLE 

A Determination of the External Forces Required 
to Move the Benchmark Active Controls Testing Model 


6. AUTHOR(S) 

Jonathan D’Cruz 



3. REPORT TYPE AND DATES COVERED 

Technical Memorandum 


5. FUNDING NUMBERS 

WU 505-63-50-15 


7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 

NASA Langley Research Center 
Hampton, VA 23681-0001 


8. PERFORMING ORGANIZATION 
REPORT NUMBER 


9. SPONSORING /MONITORING AGENCY NAME(S) AND ADDRESS(ES) 

National Aeronautics and Space Administration 
Washington, DC 20546-0001 


10. SPONSORING /MONITORING 
AGENCY REPORT NUMBER 


NASA TM-107743 


11. SUPPLEMENTARY NOTES 

This work was carried out while the author was attached to the Aeroservoelasticity Branch in the 
Structural Dynamics Division of the NASA Langley Research Center, under the auspices of 
The Technical Cooperation Program/Sub-Group 1H (Aeronautics), Technical Panel 8 (Structures 


12a. DISTRIBUTION /AVAILABILITY STATEMENT 


12b. DISTRIBUTION CODE 


Unclassified - Unlimited 
Subject Category 02 


13. ABSTRACT (Maximum 200 words) 

In view of the strong need for a well-documented set of experimental data which is suitable for the 
validation and/or calibration of modem Computational Fluid Dynamics codes, the Benchmark 
Models Program was initiated by the Structural Dynamics Division of the NASA Langley 
Research Center. One of the models in the program, the Benchmark Active Controls Testing 
Model, consists of a rigid wing of rectangular planform with a NACA 0012 profile and three 
control surfaces (a traihng-edge control surface, a lower-surface spoiler, and an upper-surface 
spoiler). The model is affixed to a flexible mount system which allows only plunging and/or 
pitching motion. Contained herein is an approximate analytical determination of the forces 
required to move this model, with its control surfaces fixed, in pure plunge and pure pitch at a 
number of test conditions. This provides a good indication of die type of actuator system required 
to generate the aerodynamic data resulting from pure plunging and pure pitching motion, in which 
much interest has been expressed. The analysis makes use of previously obtained numerical 
results. 


14. SUBJECT TERMS 

Benchmark Active Controls Testing Model, Benchmark Models Program, 

force identification, forces required to move wind-tunnel model, 

inverse problems, pure pitch, pure plunge, reconstruction problem 

17. SECURITY CLASSIFICATION I 18. SECURITY CLASSIFICATION I 19. SECURITY CLASSIFICATION 
OF REPORT OF THIS PAGE OF ABSTRACT 


Unclassified 


NSN 7540-01-280-5500 


Unclassified 


15. NUMBER OF PAGES 

30 


16. PRICE COOE 

A03 


20. LIMITATION OF ABSTRACT 


Standard Form 298 (Rev 2-89) 

Prescribed by ANSI Std 239* tB 
298-102 


















